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Abstract 



CO | We study the Lagrange formalism of the (rational) Calogero-Moser (CM) system, both 

in discrete time as well as in continuous time, as a first example of a Lagrange 1-form 
structure in the sense of the recent paper [TU] . The discrete-time model of the CM system 
was established some time ago arising as a pole-reduction of a semi-discrete version of the 
KP equation, and was shown to lead to an exactly integrable correspondence (multivalued 
map). In this paper we present the full KP solution based on the commutativity of the 
discrete-time flows in the two discrete KP variables. The compatibility of the corresponding 
Lax matrices is shown to lead directly to the relevant closure relation on the level of the 
Lagrangians. Performing successive continuum limits on both the level of the KP equation 
as well as of the CM system, we establish the proper Lagrange 1-form structure for the 
continuum case of the CM model. We use the example of the three-particle case to elucidate 
the implementation of the novel least-action principle, which was presented in |19j . for the 



' simpler case of Lagrange 1-forms. 

1 Introduction 

The Calogero-Moser (CM) model [2H] is an integrable one-dimensional many-particle system 
with long-range interactions, originally given as a continuous system with pairwise inverse 
square potential, but which also has been generalized to the trigonometric [5] and the elliptic 
case U), CI - a l so ED> an d later to a relativistic model (the Ruijsenaars-Schneider model [H[5]). 
The model has been extensively studied in both the classical and quantum cases [5] , and is 
integrable on both levels. From a physical perspective the emphasis is on the dynamics of 
real particles subject to pairwise repulsive potentials. However, in recent years the study 
of CM models has attained a wider mathematical significance, from the seminal paper [TU] 
considering the symplectic geometry associated with the model, to more recently its role in 
the representation theory of Lie and quantum algebras [11) . In this context, the investigation 
of CM systems with attractive potential, allowing particles to collide, has led to interesting 
new perspectives from the point of view of algebraic geometry, cf. e.g. |15j . 

An integrable discrete-time version of the (rational) CM model was presented in [18 , 
where it was obtained from a semi-discrete Kadomtsev-Petviashvili (KP) type equation with 
two discrete and one continuous independent variable. The construction, which we summa- 
rize in Section 2, follows closely Krichever's pole-reduction of the continuous KP equation 
leading to a connection with the CM system [13]. In the discrete case, this leads to a rather 
complicated system of ordinary difference equations (OAEs) which constitutes an integrable 
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Lagrangian correspondence, i.e., a multivalued symplectic map [16] . The construction pro- 
vides a Lax pair, whilst the classical R matrix is esentially the same as for the continuous 
case (as given by [13]), from which the Liouville intcgrability (in the sense of [TB]) follows. 
The structure of the Lax matrices was used to obtain the solution of the initial- value problem 
in an adaptation of the standard way, cf. [5] 110) . It was shown that a "naive" continuum 
limit of the discrete equations reduces them to the equations of motion of the continuous CM 
system in the case of an attractive, rather than repulsive, potential. In spite of the fact that 
this makes the discrete-time CM model perhaps less relevant from a physical perspective, it 
is certainly of interest from the mathematical point of view in that it provides us an insight 
into the basic structures underlying discrete integrability. Thus, in the same vein, we will 
use the model here to understand better the novel insights from our recent paper [19] of 
the Lagrangian multi-form structure underlying systems which are integrable in the sense of 
multidimensional consistency. 

Our focus in the present paper on Lagrangian structures has many motivations, most 
notably the possibility of quantizing the discrete-time model via a path integral formalism. 
Furthermore, whilst many integrable discrete systems admit a Lagrangian description, in 
the discrete-time case the Hamiltonian no longer provides such a natural framework as it 
does in the continuous-time case. Recently, two of the authors observed a new fundamental 
property of Lagrangians for integrable (in the sense of multidimensionally consistent) systems 
which reflects the multidimensional consistency of the lattice systems [19] . This is given by 
a closure relation which holds for Lagrangians when the system is embedded in a higher- 
dimensional lattice (meaning that the same equations hold in a multitude of sublattices of 
the multidimensional lattice together with their corresponding Lagrangians), and leads to the 
interpretation of the Lagrangians as closed forms on the extended lattice. This observation, 
based on the elaboration of various explicit examples, prompted the proposal of a new 
variational principle for integrable systems which involves explicitly the geometries in the 
space of independent variables. In [T!5] we focused on establishing and interpreting the closure 
relation of Lagrangian 2-forms for cases of integrable discrete equations with two independent 
variables, as well as of their continuous analogues, for lattice equations in the ABS list 
|20) . Subsequently, the closure property was also established for Lagrangians describing 
multicomponent lattice systems in the so-called lattice Gel'fand-Dikii hierarchy, [3T], and for 
Lagrangian 3-forms in the case of the discrete bilinear KP equation in [32], whilst in [25] 
it was shown that all equations in the ABS list possessed this property (cf. also [3T] for 
a "universal" Lagrange multi-form description for quadrilateral affine linear equations and 
their continuous counterparts). An important case which has not yet been covered is that 
of Lagrangian 1-forms, which applies to the case of integrable ODEs. Natural candidates for 
exhibiting such a structure are integrable finite-dimensional many-body systems, such as the 
Calogero-Moser systems. 

Thus, in this paper we set out to establish the Lagrange 1-form structure of the discrete- 
time CM system in the rational case. Whilst some results demonstrating that the considera- 
tions extend to the hyperbolic/trigonometric and elliptic case as well are given in Appendix 
iBl we prefer to work with this simplest case, which is most transparent, and for which we 
have exact solutions at our disposal. Starting with the discrete-time case, we will establish 
the Lagrange structure also for the continuous case by continuum limits, rather than using 
the (known) Hamiltonians for the higher CM flows. This is necessary because, as it turns 
out, the proper Lagrangians cannot be obtained by performing Legendre transformations 
separately on each of the Hamiltonians. This leads, beyond the second-order flow, to La- 
grangians containing rather complicated algebraic expressions of the higher-time-derivatives. 
Instead, we obtain a hierarchy of mixed polynomial Lagrangians in terms of these higher- 
order derivatives from the discrete-time case by systematic expansions (indicating en passant 
that Legendre transformations should probably not be implemented order-by-order in the 
flows but rather altogether in one stroke). In all these considerations, the link with the un- 
derlying KP system is instrumental in deciding how to perform the higher-order continuum 
limits by systematic expansions and in terms of the proper parameters. 

The organization of the present paper is as follows. In Section 2, we present the full solu- 
tion of the semi-discrete KP equation, which requires commuting flows in different discrete- 
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time directions, and thus generalizes the one-dimensional discrete-time CM flow in [18]. In 
Section 3, we establish the closure relation of the discrete-time Lagrangian 1-form of the 
discrete-time CM model through the compatibility of the time-part Lax matrices, which are 
of Cauchy type. In Section 4, we investigate the semi-continuum limit, or skew limit, of the 
semi-discrete KP equation. This limit leads to an equation which we refer to as the semi- 
continuous KP equation (defined in terms of two continuous and one discrete independent 
variable), from which we obtain by pole-expansion the corresponding differential-difference 
system which we coin the semi-continuous CM system. In Section 5, the latter, which acts 
as a generating system for the CM hierarchy, is a Lagrangian system in its own right obeying 
a closure relation of differential-difference type with the original discrete-time Lagrangian of 
the CM system. In Section 6, the full continuum limit is performed, and we thus recover the 
usual fully continuous CM hierarchy in mixed form, together with continuous Lagrangians 
exhibiting a closed 1-form structure in Section 7. We concentrate in Section 8 on the 3- 
particle case, governed by the CM flow and its first higher-order counterpart, in order to 
implement the geometric variational principle which was formulated in [15]. Thus, we show 
how the principle leads in this case to the derivation both of the closure relation as well as 
the pertinent Euler-Lagrange (EL) equations. We end the paper with some discussion of 
open problems and possible extensions. 



2 Pole reduction of the semi-discrete KP equation 

In this section we review the connection between the semi-discrete KP equation, following 
on from [12] [17], and the discrete-time CM model as in [18, 24 . This gives us an occasion 
to introduce appropriate notation which we will use throughout the paper. Furthermore, it 
helps us to identify the commuting flows needed to establish the corresponding nontrivial 
solution of the semi-discrete KP equation. 

The semi-discrete KP equation that was used in |18] is 

d^(u — u) = (p — q + u — u)(u + u — u — u) , (2-1) 

where p and q are two lattice parameters, u is the classical field, and u represents the discrete 
shift of u corresponding to a transition in the "time" direction, while u represents the discrete 
shift of u corresponding to a translation in the "spatial" direction. There are two discrete 
variables n, m to which the lattice parameters p, q and shifts u, u correspond respectively, 
and one continuous variable £. 

We observe that eq. (|2.1j) is a consequence of the compatibility condition of a scalar Lax 
pair which reads 

4> = <t>i + (P + u - u)(j) , (2.2a) 
4> = (j)c + (q + u - u)<p . (2.2b) 

Taking u to have the form 

N 1 

u = V- 7 -, (2.3) 

—-[£,- Xi{n,m) 

we find that the corresponding solution of the linear Lax equations is 

*-('-s|{^o)c+*)"(«+*r« ,< . <"> 

where k is a new spectral parameter and the bi are yet to be determined. Inserting expressions 
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(|2.3[) and (12 .4[) for u and 4> into the first Lax equation (|2.2a|l . we obtain the relations 



* +k * = HE— -E— h-E^ (2 - 5) 
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1,2,..., AT. (2.6) 



N P N V + P 

l = y -^v - y , (2.7) 

■ - , . . , 

z,j — 1 J i,3 = 1 

where £Iy are the generators of GLyv, i.e., matrices with entries (Eij)ki = SikSji, we can now 
rewrite equations (|2.7[) and (|2.8p as 

(p + k)b = ke + Lb, (2.9a) 
(p + fc)b = fee + Mb , (2.9b) 

which forms a N x 2V matricial Lax pair. The compatibility of equations (|2.9aD and (I2.9bp 
gives us the relation 

LM-ML)b + k[L-M)e 



N 

E 



* * \ ~ ~ — Tj / * * 1". — 



1 = 1 \Xi — Xl X i X l / r =1 * T 2 ^ 



k.e-Ab|=ll 

AT x, — a;, 



and since this must holds for arbitrary value of k we deduce the discrete-time equations of 
motion of an A-particle system: 



N , x N 



Vf + |-2y =0, where i = 1,...,JV , (2.10) 



in which the under-shift Xj is short-hand for a shift in the variable n over one unit in the 
nagative direction, i.e., Xi(n,m) — Xi{n — l,m). Furthermore, if (|2.10p holds, both the 
following matrix relations are satisfied: 

LM = ML , and (L - M)e = , (2.11) 



and consequently the inhomogeneous Lax system (|2.9[) leads to an isospectral problem in 
terms of the Lax matrices L and M . As was shown in [IS] eq. (|2. 10[) is an integrable 
symplectic correspondence, in the sense of Liouville, |16j . whose exact solutions Xi(n,m), 
which denote the position of the particles Xi at the nth time step (with fixed to) , are obtained 
by a linearization procedure implementing the Lax pair (for details, see Appendix |A~| . The 
resulting solution of the initial value problem, imposing initial values Xi(0, to) and X{{1, to) = 
Xi(0, to), can be obtained by solving the secular problem given by the characteristic equation 
of the matrix 

Y(n, m) = Al n Y(0, m)A£ - nA^ 1 . (2.12) 
subject to the constraint on the initial value matrix 

[Y(0,m) , A L ] = 7 + rank 1 , 
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(in which I denotes the N x N unit matrix). The eigenvalues of Y(n, m) given by (|2. 12|) are 
the particle positions Xi(n,m), and they are determined up to permutations of the particles, 
which accounts for the multivaluedness of the system of OAEs. Here the A^ denotes the 
diagonal matrix of eigenvalues of the Lax matrix L and the matrix Y (n, m) is related to the 
the diagonal matrix X(n,m) — 'Y^ =l Xi{n,m)Ea of the particle positions by a similarity 
transformation. It follows from (|2.11l) that the matrix Al is invariant under time shifts, i.e., 
A L = A L . 

We note that (|2.12[) resolves only the n-dependence of the solution of the KP equation, 
and to provide the full solution we need also to consider the dependence on m. Obviously, 
the latter should come from the consideration of (|2.2b[) . which can be treated in much the 
same way, and lead to the following relations: 

(q + k)b = ke + Kb, (2.13a) 
(q + k)b = ke + Nb, (2.13b) 

where K and N take exactly the same form as L and M, with just the "~ " shift replaced 
by the """* " shift. The compatibility condition leads to 

KN = NK , and (K - N)e = , (2.14) 

and the equations of the motion can be directly obtained from (|2.10p by replacing the " ~ " 
shift by the "~" shift, 

N , v N I 

-2V =0, where i = l,...,N, (2.15) 

X-r X ■> 



£ 



in which the under-shift is short-hand for a shift in the variable m over one unit in the 
negative direction, i.e., Xi(n, m) = Xi(n,m — 1) The exact solution of the system of OAEs 
(I2.15[) can be obtained in the same way as (I2.12j) . namely by considering 

Y(n, m) = A K m Y(n, 0)A^ - mA^ 1 . (2.16) 

subject to the constraint on the initial value matrix 

[Y(n,0), A K ] = J + rankl , 

where the matrix Ak is invariant under time shift, so that = A^-. 

The n- and m-parts of the solutions of course can only be combined if the corresponding 
flows are compatible. In particular this requires that the Lax matrices L and K commute, 
which in turn suggests that these matrices can be simultaneously diagonalized. Further- 
more, it requires the compatibility of the flows in the two lattice directions, which leads to 
the relations 

(p-q)b = (L — K)b , (2.17) 
{p-q)ke = (M -N)ke + (M N -N M)b . (2.18) 

Eq. (|2.17[) leads to the relation 

JV 

p-q 



N ( 1 1 \ 

£ -J, i = l,...,N, (2.19) 



1 = 1 



Xi X[ X<i X\ 



whereas (I2.18P yields subsequently 



which yields in addition 



p-<i=J2 — --^^ ' i = 1 '---' iV - ( 2 - 2 °) 

j—l ™ *£j J 

Furthermore, the relation 




i,j = l,...,N, (2.21) 



which is a consequence of (|2.19[) and (|2.20p . guarantees that the zero-curvature condition 
MN = NM holds, which in turn implies that (M - N)e = (p - q)e . 

Equations (|2.19[) and (|2.20l) , which we will refer to as the constraint equations, guarantee 
that the discrete flows in the variables n and m commute, and hence that the corresponding 
linear equations can be simultaneously solved. In fact, the downward shift of (|2.20p yields 




which implies 




(2.22) 



which is also a consequence of (|2.10p and (|2.15[) . and which expresses this compatibility with 
the sets of OAEs. It is, thus, the consistency of the equations (12.12)) and (|2.16l) . that lead 
us to the complete solution of the semi-discrete KP equation. In particular, from eq. (|2.19l) . 
it follows that the diagonal parts of the matrices L and K differ by (p — q)I, and hence we 
can identify 

A L = pi + A , A K = ql + A , 

with common diagonal matrix A. Using this identification we can now combine the solution 
given by (|2.12[) for the CM system of OAEs in the n-direction and the solution (|2 . 1 6[) in the 
m-direction into a simultaneous solution which satisfies also the commutativity constraints 
(I2.19j) and (|2.20p . leading to the following statement: 

Proposition: The eigenvalues X\{n, m), . . . , xn{n 1 m) of the N x N matrix 

Y(n,m) = (pI + A)- n (qI + A)- ,n Y(0,0)(pI + A) n (qI + A) m -n(pI + A)- 1 -m(qI + A)- 1 

(2.23a) 

in which the initial value matrix 1^(0,0) is subject to the condition 

[Y(0, 0) , A] = I + rank 1 , (2.23b) 

obey both the discrete-time Calogero-Moser systems given by eqs. (|2.10p and (|2.15p as well 
as the systems of constraint equations given by (|2.19p and (|2.20p . 



The proof is presented in Appendix [X] In order to make a connection with an initial value 
problem, we mention that the initial value matrix 1^(0, 0) can be obtained from the diagonal 
matrix of initial values X(0, 0) by a similarity transformation with a matrix [7(0,0) which 
is an invertible matrix diagonalizing the initial Lax matrices L(0, 0) and K(0,0). To find 
the latter, we need the initial values 2:^(0,0), Xi(l,0) and Xi(0, 1), (i = 1, . . . , N). We note 
that the secular problem can, hence, be reformulated as one for the following matrix 

X(n, m) = X(0, 0) - nL _1 (0, 0) - mK-\0, 0) , (2.24) 
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and hence the solution is provided by the roots of the characteristic equation: 

N 

Px(x) = det(xJ — 3C(n, m)) = Y\_( x ~ x i( n > m )) ■ (2.25) 

i=l 

Returning now to the solution of the semi-discrete KP equation (|2.1I) , inserting the roots 
obtained from (|2.25[) into the expression (|2.3[) we obtain the required pole solutions. Fur- 
thermore, for the solution of the corresponding eigenfunction </>(n, m) as given in (|2.4[) for the 
eigenvalue problem, the vector b(n, m) must be determined. This is obtained, for arbitrary 
value of the spectral parameter k, from (|2.9a[) . or equivalently (|2.13a|) . and is given by 

b = (p + k — L) 1 ke , or equivalently b = (q + k — K) 1 fee . (2.26) 

Plugging all the results into (|2.4I) . we obtain the explicit form of the eigenfunction <f)(n,m) 
solving eqs. (|2.2a|) and (|2.2bl) . Furthermore, eliminating the function u from the latter, we 
obtain the following nonlinear differential-difference equation for cf> itself: 



<t>(cl><f> e -<f><P e ) = (<t><P-(t>(t>)((t>-<t>) , (2.27) 
whose pole-solutions are given by the construction outlined above. 



3 The Lagrangian 1-form and its closure relation 

In this section, in order to derive the Lagrangian 1-form structure, we first establish a direct 
connection between the Lagrangian of the discrete-time CM model, as was presented in [T5] . 
and the temporal part of the Lax representation. We will show that this connection, through 
the compatibility of the matrices M and N of the previous section, will immediately lead 
to a closure relation for the corresponding Lagrangians. 

Let x = (#1, X2, ■ • • , xn)- Equation (12. 10)) can be computed from the variation of a 
discrete action given in [THj as 



S(„) [x(n, m)] = y^£( n )(x,x) = ^ - log \xj - Xj \ + ^ log la;. 



N N 

X, 



\ 



(3-1) 



where £ir n \(x, x) is the Lagrangian corresponding to the "~ " direction, and the sum over n 
represents the sum over all discrete-time " ~ " iterates. Similarly, the action corresponding 
to the " ^ " direction takes the form 



N N 



S(m) [aj(n, m)] = y^£( m )(x,x) = | - \og\xt -Xj \ + log \xi - Xj \ | , (3.2) 

where Zi m \{x,x) is the Lagrangian corresponding to the "~ " direction, and the sum over 
m represents the sum over all discrete-time " iterates. 

The matrices M and N are Cauchy-type matrices, and so the exact forms of the deter- 
minant can be written as: 
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and we also have that 
log|det(M)| 



iV 



A' 



= ^2 (log \xi -Xj\ + log \xi - Xj\) - ^ lo S \ x i ~ x i\ > (3.5a) 



'■j=i 

N 



N 



log |det(AT)| = 



^ (log \x{ -Xj\ + log \xi - Xj I) - ^ lo S l x i _ %l ■ (3.5b) 



'■,i=i 

i<3 



Remarkably, the action (|3.1|l can be obtained by considering the infinite chain product of 
the matrix M in the following way 



>(n) 



log 



det Yl M ( n ) 



Vn— +oo 



A' 



E X! ( log 1^ ~ x i \ + log 1^ ~~ ^'D ~ E log l Xi ~ ^'1 



n \ i j — 1 

El - E 



A r 



+ ^ log 



\ 

/ 



E £ 

(x,x). (3.6) 



The same procedure can performed on the matrix TV in order to get the action 



>(m) 



log 



E 



det Yl Ar ( m ) 



\77l — + 00 



A 



\ 



^ (log |x< - Xj| + log \xi - Xj\) - 2J log - %l 



N 



A 



E - E log \ Xi - %i + E lo § 



E^ 



(m) I, 3 '? ' 



(3.7) 



KJ=1 



We note that the connection between the time-part Lax matrix and the Lagrangian also ex- 
ists in the cases of the elliptic and trigonometric discrete-time CM systems (see Appendix [Bj. 

We now consider the compatibility relation MN = NM, which we recall is satisfied if 
both (|2.19p and (|2.20l) hold, and we rewrite this as 

log |det(M)| + log |det(JV)| - log jdet(iV)| - log |det(M)| = , (3.8) 



or 



L {n) (x,x) - L {n) (x,x) -L (m) (x,x) +L {m) (x,x) = , (3.9) 
where we define S = x « to be the centre of mass and 

L {n) (x,x) = log|det(Af)|+p(H-S), (3.10a) 
£ {m) (x,x) = log|det(JV)|+ 9 (S-H). (3.10b) 

The last terms in (j3~T0a| and (j3.10b[) are the total derivative terms involving centre of mass 
motion which can be separated from the relative motion. The extra equation comes directly 
from ([331) 

3+5=5+1, (3.11) 
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which holds on the solutions. 



We would now like to interpret p. 91) as the closure relation for a Lagrangian 1-form. 
Let &i represent the unit vector in the lattice direction labeled by i and let any position in 
the lattice be identified by the vector n, so that an elementary shift in the lattice can be 
created by the operation n i— > n + e$ . Since the Lagrangian depends on x and its elementary 
shift in one discrete direction, it can be associated with an oriented vector e, on a curve 
r,*(n) = (n, n + e,-), and we can treat these Lagrangians as defining a discrete 1-form Li{n) 

Li(n) = £>i(x(n),x(n + e»)). (3-12) 

Proposition: As functions on the two-dimensional lattice, the discrete-time Lagrangians 
given in eq. (|3.10p satisfy the following relation 

Li(x(n + ej), x(n + + ej)) — Li(x(n), x(n + e,;)) 

—Lj(x(n + ei),x(n + ej + e«)) + £j(x(n), x(n + ej)) = . (3.13) 



Equation (|3.9p represents the closure relation of the Lagrangian 1-form. The point is 
that the closure relation (|3.9I) is derived from the Lax equation (the compatibility of the 
matrices M and TV), together with the property of determinants of Cauchy matrices. In 
contrast, in those cases of Lagrangian 2-forms and 3- forms in [1151 HIJ 122 > the equation of the 
motion must be invoked in order to show that the closure relation holds. Also in those cases 
the Euler-Lagrange equation stemming from the variational principle was not exactly the 
equation of motion needed to verify the closure relation; rather it was a discrete derivative, 
or a sum of copies, of the original equation of motion. 

Choosing a discrete curve T consisting of connected elements T^, we can define an action 
on the curve by summing up the contributions from each of the oriented links Ti in the 
curve, to get 

S(a!(n);r) = ^2jC i (x{n),x{n + e i )). (3.14) 

The closure relation p.9[) is actually equivalent to the invariance of the action under local 
deformations of the curve. To see this, suppose we have an action S evaluated on a curve 
r, and we deform this (keeping end points fixed) to get a curve T' on which an action §' is 
evaluated, such as in Figure [TJ 



Figure 1: Deformation of the discrete curve T. 

Then S' is related to § by the following: 

S' = § — Li(x(n + ej),x(n + ei + ej)) + Li(x(n),x(n + ei)) 

+Lj(x{n + ei),x(n + ej + e^)) — Lj(x(n),x(n + ej)). (3.15) 

9 



Equation (|3.15[) shows that the independence of the action under such a deformation is locally 
equivalent to the closure relation. The invariance of the action under the local deformation 
is a crucial aspect of the underlying variational principle. 

m m 




Figure 2: The effect of changing variables on the discrete curve I. 



Next, it is interesting to investigate further how to derive the discrete Euler-Lagrange 
equation from the variational principle. The system given in Fig. (TTJ) is difficult to work with 
because of a matter of notation. For general discrete curves it is cumbersome to implement 
the variational principle because of the notation it would require. We will, however, demon- 
strate how the principle works for a few simple cases: (a) the curve shown in Fig. ([2]), (b) 
the curve shown in Fig. 

(a): The curve shown in Fig. We now introduce a new variable N = n + m (which 

will play an important role in the next section) together with the change of notation 

x(n, m) h-» x(N, m), x := x(N + 1, m) and x := x(H, m + 1) , 

and so we work with the curve given in Fig. (JSfr)- The action evaluated on this curve can 
be written in the form 

mi— 1 'mi — 1 

S[x;r']= -£ (N )(a:(No 1 m),s(No-l,m))+ £ L {m) (x(N Q -l,m),x(N ,m + l)p.W) 

m—mo m—mo 

where 

N (\ 1 \ N 

£(N) 0, y) = ( 2 log ~ y i I + 2 log ' Xi ~ Xj ' J ~ log l yi ~ x i I 

N 

(3-17) 

i=l 

N fl 1 \ N 



2 

AT 



+p£><-w). (3.18) 

i=l 

The minus sign in p,16[) indicates the reverse direction of the Lagrangian £(m along the 
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horizontal links. Performing the variation x t-¥ x + Sx, we have 



5S = = 

"x^ 1 f 9£(N)(x(N ,m),£c(N -l,m) 
2^ I a.^ M „^ <5£c(N ,m) 

m—mQ ^ 



9a;(No, m) 



d£( N) (a;(No,m),3;(No - 1, 
9a;(No — 1, m) 

/&C (m) (a;(N - l,m),x(No,m - 1) 
- 2^ I HZTm — ~ I T\ dx{N ,m + l) 

m-mo ^ 



<5a;(No — 1, m) 



dx(N ,m + 1) 

<9£ (m) (a;(No - l,m),x(N ,m + 1) 



+ - 



9a;(No — 1, m) 

We now obtain the Euler-Lagrange equations 

9£ (N )(a;(No,rn),a;(No - l,m) <9£ (m )(a;(No - l,m - l),a;(No,m) 

9a;(No,m) 9a;(No,m + l) 

5£( N) (a;(No,TO),a;(No - l,m) S£ (m )(x(N - l,m),x(N ,m+ 1) 



&r(N - l,m) (3.19) 



which produce 



p-q 



p-q = 



dx(N„ - l,m) 



N 

E 

N 

E 



&r(N - l,m) 



= 0, 
= , 



1 



1 



Xi(U ,m) - Xi(H + l,m - 1) Xj(N ,m) - Xj(N - l,m) 
1 1 



V^iC^o,™) - a=i(N + l,m) x,(No,m) - aJ 4 (N - l,m + l) 
which are equivalent to (|2.19p and (|2.20[) . respectively. 



(3.20a) 
(3.20b) 

(3.21a) 
(3.21b) 



mi 




m 



(n ,m ) "a 
(a) 



mi 



(N' ,m ) N'^1% + 1 
(b) 



N' 



Figure 3: The effect of changing variables on the discrete curve II. 

(b): The curve shown in Fig. (|3]): Introducing the variable N' = n — m, the corresponding 
curve is given in Fig. (JSJj). The action evaluated on the curve V reads 



mi — 1 



mi — 1 



§[x;T'}= £ ( N')(*(No,ro),aj(N{, + l,m))+ £ L (m) (x(N' + l,m),x(N' ,m+ 1)#.22) 
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where 



N (\ 1 \ N 



N 



+p^2(x i -y i ), (3.23) 



i=l 



£(m)(a:,y) = ( ^ lo s\xi-x j \ + -log\y i -y j \\ - ^ lo. 



N 

+pJ2(xi-Vi)- (3-24) 



Performing the variation x t-^- x + 5x, we have 
5§ = = 

S I Rn^) ( °' m) 

5£ (N , ) KN»,x(Nj + l, m ) 
+ ^(N'o + l.m) ^(N + l,m, 

+ £ { dx(U> ,m+l) MNo,m + l) 

m=m x v u ' 

aC (m) (g(N{, + l,m),s(N{,,m + l) , 
+ to(N^ + l, m ) **(No + l,m)]..^:„ 

We now obtain the Euler-Lagrange equations 

d£ (NI) {x(N' ,m),x(K + l,™) d£ {m) (x(N' Q + l,m-l),x(K,m) 

dx(N' ,m) dx{N' ,m + l) 

dL m (x(N' ,m),x(% + l,m) dL (m) (x{N' + l,m),x(N' ,m + 1) 



, (3.26a) 
= , (3.26b) 



dx(N' + l,m) dx(N' + l,m) 

which produce 

N / I i 

q ~ P = ^ Ui(N(,,m) - x~(% - l,m+ 1) " ^(N^m) - ^(1% + l,r^ ) ' (3-27a) 



1 1 



P-9 = 



^ (ai(N{, + l,m) - b<(N{, + 2, m) x<(N{, + l,m) - Xi(N' , m + 1) / ^ 



which are again equivalent to (|2.19D and (|2.20[) , respectively. 

By working with the specific type of curves given in Fig. ([2]) and Fig. <j3j> , we can perform 
the variational principle with either the new variables (N,m) or (N',m). Even though we 
have these two possibilities, the natural choice is to work with (N,m). This is because the 
variable N' leads to difficulties in performing the skew limit in the next section. This can, in 
fact, be seen from the structure of the plane wave solution (|2.4[) . 

It is not yet clear whether the processes of varying the curve of the independent variables, 
from which we get the closure relation, and that of varying the dependent variables, from 
which we get Euler-Lagrange equations, commute. We defer consideration of this point to a 
future publication. 
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4 The semi-continuum limit: skew limit 



In this section, we study a continuum analogue of a previous construction in Section 2 by 
considering a particular semi-continuous limit. Since the semi-discrete KP equation (|2.1[) 
contains two discrete variables n and to, we could perform a continuum limit on one of 
these variables separately, while leaving the other discrete variable intact, and thus obtain 
a semi-continuous equation with one remaining discrete and two continuous independent 
variables. Alternatively, we can first perform a change of independent variables on the 
lattice and subsequently perform the limit on one of the new variables. The advantage of 
the latter approach over the former is that it often leads in a more direct way to a hierarchy 
of higher order hows. Adopting the latter approach in this section, we use a new discrete 
variable N := n + to, and perform the transformation on the dependent variables by setting 
u(n,m) i — y u(N,to) =: u, which leads to the following expressions for the shifted variables: 

u = u(n + l,m) i — ^ u(N + 1, m) =: u , 
u = u(n,m + l) i — ^ u(N + 1, m + 1) =: u , 

u = u(n+l,m + l) h» u(N + 2, m + 1) =: Tj . 

Rearranging the discrete exponential factor in (|2.4[) , we have 

s ( 1 bj{n,m) \ N / q + k \ m k , 

0(n,TO = I-7V7 7 T )(P + k ) — TT e > 

V fci^ - Xi(n,m) J \p + kj 

We perform the limit n — > — cxd, to — > 00, e — > while keeping N fixed and setting e = p — q, 
such that £to = r remains finite. Focusing on the penultimate factor in (|4.ip we have that 

lim (1 = lim (1 T — ) =e - ^, (4.2) 

m-»oo y p + k J ?ti->oc y to(ji + k) J 

so that the solution to the Lax equations (I2.2a[) and (|2.2b|) takes the form 

wr K , ~*t«^) < ' + * ,v<e ~*- <t3) 

We would now like to see the effect of this limit on the semi-discrete KP equation. Wc 
rewrite eq. (12.11) as 

ln(p — q + u — u) =u — u — u + u. (4-4) 

Taking the limit, we have 

<9{ ln(e + u — u) = u — u — u + u. (4.5) 

Setting r = tq + me, where eo is a background constant of the continuous variable, and 
applying the Taylor expansion, we get 

u(r + e) = u(r) + £^u(t) + i^^ufr) + . . . (4.6) 

where u = ^u. The leading order 0(5°) gives us 

<%ln(l + u) = u+u-2u. (4.7) 

Eq. (|4.7p . which contains two continuous and one discrete variable, we call the semi- 
continuous KP equation (to distinguish it from the semi-discrete KP equation, which contains 
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one continuous and two discrete variables). 

Remark: The connection to the Toda system. 

Introducing the variable 1 + u =: cxp(v) wc can write eq. (|4.7I) in the form 

v T £ = 2cxp(v) — cxp(y) — cxp(v). (4-8) 

Introducing next another variable y through the relation v = y — y, wc obtain 

Yr? - Yrt; = (exp(y - y) - exp(y - y)) + (cxp(y - y) - exp(y - y)) , (4.9) 

from which we can identify 

y r? = cxp(y - y) - cxp(y - y) . (4.10) 

This is the well-known two-dimensional Toda lattice, cf. [25]. A Lagrangian exists for this 
system is given by 

£>2D-Toda = T^y? + eY _ y (^.11) 

To follow the reductions of the Lagrangian structure for the CM system, it would be desirable 
to have a Lagrangian for the semi-discrete KP equation (I2.1j) and for the semi-continuous KP 
(|4.7[) equation, but these are strangely elusive. Obviously, the fully continuous KP equation 
does possess a Lagrangian structure which is easy to establish by inspection, and furthermore 
the 2D Toda lattice (|4.10l) reduces to the KP equation through some specific continuum lim- 
its, see e.g. [23]. Thus it is conceivable that CM lagrangians could be established through 
the "Toda route", but we will not pursue this line in the present paper. We would like to 
mention also Krichever's elliptic analogue of the one-dimensional Toda lattice [32], which 
remarkably contains a term corresponding to the elliptic discrete-time CM given in [24] . We 
will give some results on the Lagrangian 1-form structure of the discrete-time elliptic CM 
model in Appendix [B] but our main concern in the present paper is to study the Lagrangian 
1-form structure in the simplest possible case, for the sake of transparency. 

The skew limit on the Lax equations: To obtain the Lax representation for the semi- 
continuous KP equation (|4.7p . we perform a similar limit on the Lax equations (|2.2ap . (|2.2b[) 
by making the transformation cf>(n,m,£,) H > 0(N,m,£) with similar replacements as above 
for its lattice shifts, i.e. 

<f)(n+ l,m,£) I-+ (f>(N + l,m,£) =: 4> , (j){n,m + 1,£) ■->• <f>(N + l,m + 1,0 =: 4> . 
Thus, we obtain 

4> = c + (p + u-u)& (4.12) 
4> = <fe + (g+u-u>. (4.13) 
Using a Taylor expansion as in the previous case, we can write (|4.13[) as 

4> + e^+... = S + (>-£ + u-u-eu' + ...)<£, (4.14) 
where = -§^(j>- Then to leading order 0(e) we have 

0=-(l + u)0. (4.15) 

Equation (|4. 1 5[) is a mixed equation, with one discrete and one continuous variable. It is 
easy to check that eq. (|4.7I) arises as the compatibility condition of the Lax pair consisting 
of gUP and (J132J. 
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Inserting the special form of the solution, eq. (14.3[) . into eq. (14.12|) . we recover the set of 
equations f|2.5[) and (|2.6I) . obviously with the replacements of Xj(n, m) by Xi(N,m) and with 
Xi = Xi(N + l,m). For the sake of selfcontainedness we write the Lax corresponding matrices 
(which are essentially the same as before) in the new notation as follows: 

N F N F 4- F 

— 1 J 1,J— 1 

Ej, 



The compatibility between the matrices L and M produces the equations of motion 

N i 

-2^ =0, where i = l,...,iV, (4.18) 




3=1 
jjii 

implying as before the system of equations (|2.15[) . but now viewed as part of a semi- 
continuous Calogero-Moser system. This set of equations will be complemented by equations 
involving derivatives with respect to the variable r, which are obtained by inserting the so- 
lution (|4.3p into eq. (|4.15[) and obtaining the complementary equations for b^. Equating to 
zero the resulting coefficients of (£ — x^) -1 , (£ — x^)~ 2 and (£ — x$) , we obtain 



N 



{p + k)b % 



^(x,-x,-) 2 ' 



(p + k)bi = 



3 = 1 



N 

E 

3 = 1 



(4.19) 
(4.20) 

(4.21) 



Introducing the matrix 



N 

E 



- B- ■ 

2^3' 



(4.22) 



together with the matrix M as given above, we can rewrite equations (|4.19p and (|4.20l) as 

(p + k)b = b + Ab, (4.23) 

(p + k)b = ke + Mb. (4.24) 

Note that (|4.23[) and (|4.24p can be directly obtained by taking the skew limit on (|2.13bl) in 
order 0(8°) and 0(S) respectively. The compatibility condition between (|4.23[) and (|4.27al) 
leads to 



(p + fc)(A-M-MA M b + k (I + M A M ) e = 



which, using (|4.22f) and (|4.17p . leads to 

t N 

E 



N 

E 

i,3=l 



Ei 



■I'l 



E 



M _1 ((p + fc)b-fce) 



N 



X, 



^ (x, - x ( ) : 



Eye 
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This, using (|4.21l) . implies that we also have 



N 



j=l *"3> 

and thus, combining these two relations, we get: 




(Xi-Xj) 2 (Xj-Xl) 2 



(4.26) 



We observe that if we use (|4.25p . we recover (|4.21[) from (|4.26[) . Furthermore, we can show 
that eq. (|4.26l) is a consequence of the skew limit of the eq. (|2.21[) of order 0(8). On the 
level of the matrix relations, this implies 

(l + A-M)e = 0, (4.27a) 
MAM 1 A + M = 0. (4.27b) 

Remark: We observe that (|4.21l) and (|4.25[) can be obtained by taking the skew limit on 
(12"H?|) and ([2T2U| in order 0(5), respectively. 

The skew limit on the solution: we finish this section by performing the skew limit 
on the full solution (|2.16[) . Let us observe that due to the fact that the matrices L and K 
are related linearly through the constraint (12.191) (this relation amounts to a shift over (p — q) 
times the unit matrix) we can identify the diagonal matrices of eigenvalues as follows: 

A L =p + A, and A K = q + A , (4.28) 

where A is the matrix which is independent of both p and q (i.e., independent of the direction 
of the lattice). Then the full solution (|2.16p can be expressed in the form 

Y(n,m) = (p + A)- n (q + A)- m Y(0,0)(q + A) m (p + A) n 

-n(p + A)- 1 - m(q + A)- 1 . (4.29) 

Using the definitions of the variable N = n + m and p — q = e, we obtain 

Y(N,m) = (p + A)- N (l-e(p + A)- 1 y m Y(0,0)(l-e( P + A)- 1 ) m (p + Af 

-N(p + A)- 1 + m ((p + A)- 1 - (q + A)- 1 ) . (4.30) 

Taking the limit, we have 

lim Y(N, m) = Y(N, t) = (p + A)- N e r(p+Arl Y(0, 0) e - r(p+Arl (p + A) N 

m— >oo 

-NCp + A)- 1 -t( P + A)- 2 . (4.31) 

This equation represents the full solution after taking the skew limit. The positions of the 
particles Xi(N,m) can be determined by computing the eigenvalues of (|4.3ip . 



5 Semi-continuous Lagrangian and closure relation 

Having obtained, in the previous section, the skew limit of the Lax equation and the corre- 
sponding differential-difference system comprising eqs. (|4.21l) and (|4.25[) . together with eq. 
(I4.18p . we now proceed to present the corresponding Lagrange form. First, we observe that 
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eq. (I4.18|) can be once again be obtained by implementing the usual variational principle on 
the following action S(n) given by 

( N N (\ 1 

§(N) [x(N, t)\ = £n = [ - zl log ' Xi ~ X J I + H ( 2 l0g ~~ Xj ' + 2 l0g ' Xi ~ Xj 

N N \ < 'J =1 *^ 

AT \ 

+p^(x l -X J ) (5.1) 
i=l / 

where now the Lagrangian £(n) involves variables X; shifted in the discrete variable N instead 
of the original variable n, and the corresponding discrete Euler-Lagrange equation reads: 

yielding (|135|) . 

Next we observe that the combined differential-difference equation (|4.26[) can be obtained 
from the two-dimensional action 



S (t) [x(N,t)]=]T / dr 

M « 



x 4 



Xj — X 



yielding fl4.26|) from the two-dimensional EL equations 

dxj dxi dr \dxi J ' 

However, we will argue that this is not the correct point of view to take in the context of 
the CM systems, because they consist after all of ODE's rather than PDE's or differential- 
difference equations. In fact, this action, containing a summation over the variable N and 
an integration over the variable r can never be obtained as a continuum limit of the one- 
dimensional discrete actions presented in Section 3, which contain one single summation. 



The skew limit on the action: We note that the actions in (|5.1[) and (15. 3[) arc obtained 
directly from the corresponding equations of motion (|4.18p . (I4.21[) and (|4.25l) . respectively. 
An interesting question is whether we can obtain the mixed discrete and continuous 1-form 
structure directly by taking the skew limit on the fully discrete 1-form action (|3.14p . This 
will be rather tricky because it will sensitively depend on parametrization of the curve. Fur- 
thermore, in order to perform a systematic expansion on the action functional one would 
need to employ not only a Taylor expansion on the relevant lattice shifts, but at the same 
time also the opposite of the Taylor expansion on the sums. For the latter we can employ 
the following "anti- Taylor" expansion of the sum of a discrete function f(n) =; F(t), where 
t = To + ne with To fixed: the formula 

£ f(n) - \ f 2 dr F(r) - \ (F(r 2 ) - F(n)) + ± ((F'(r 2 ) - F'(n))) + . . . , 

n— ni Tl 

which guarantees that the analogue of the fundamental theorem of integration holds, namely: 

Tl2 — l 

E (f( n + !) - /(»)) = /M ~ f(ni) - F(t 2 ) - F( n ) . 

n—n\ 

We now would like to pursue this task by considering a simple discrete curve Y given in Fig. 
The action corresponding to the curve has been given by (|3.14[) . Performing the skew 
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limit, the curve has been transformed to Fig. (JJJd) with the new variables (N, m). Using the 
anti- Taylor expansion, the action (|3.16[) becomes 

S[x(N,m) ; r'] = ~ J "rfr£ (m) (x(N ~l,r),x(N ,r-£))-i[£ (m) (x(N -l,T 2 ),x(No,T 2 - £ )) 

- £ (m) (x(N - l,Ti),x(N 0) Ti - e))] + 

\J 2 dr£ (N) (x(N ,r),x(N -l,r)) + i[£ (N) (x(N ,r 2 -e),x(N -l,T 2 )) 

- £ (N) (x(N ,T 1 -£),x(N -l,r 1 ))] + (5.5) 

We also use the fact that the Lagrangian is an antisymmetric function under interchange of 
the arguments 

£(x(No,r) ) x(No + l,r)) = -£(x(No + l,r) ) x(N ,r)) , 
together with expansion with respect to e, then the action becomes 

5x(N 0) r)' 



§[x(N, m); r'] = J 2 drL (T) ^x(N - 1, r), ■ 



(5.6) 



where &{ T ) is given in 

By using a specific choice of curve given in Fig. ([2]), we can perform the skew limit on 
the discrete action leading to what is called the "semi-continuous" action in 



The skew limit on the discrete-time closure relation: We would like to consider 
the closure relation of the semi-continuous Lagrangians. Taking, as we did in the equations 
of motion, the skew limit on the discrete-time closure relation eq. p.9[) . implementing a 
change of variables on the lattice first, we obtain: 

£(N) — £(N) = <C(m) ~ <£(m) => 
(<£(N) + e£(N) + —J - <£(N) = (^(N) — £ ^(t) + — (^(N) — e ^(t)) ■ 

Equating the coefficients of e gives us 



d£(N) 
dr 



£(r) ~ £ (t) , (5.7) 



where the Lagrangian £( T ) is obtained by performing the skew limit on the Lagrangian £>( m ) 
in eq. (|3.10bl) to give 



N ~ N ~ 

X-i" V ^ X )' 



^)=Tr(M-A)- S + § = ^^- + ^-^-- S + g-pg ) (5.8) 

which produces the semi-continuous equations of motion. The derivation of (|5.8p is given in 
Appendix [Cj 

Furthermore, the variational derivatives with respect to x(t) and x(r) yield: 

= and — ^ = , 5.9 

which yield the equations 

dL {T) 



N 



dL [T) 



3- 
N 
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which are identical to eqs. (|4.21[) and (|4.25|) . 



The relation (|5 . 7[) represents a closure relation between the semi-continuous Lagrangian 
£ii t ) an d the semi-discrete Lagrangian £(n). This points again to a Lagrangian 1-form struc- 
ture but now of a mixed, part discrete and part continuous, type where the total action is a 
combination of a discrete part of the action S(n) and continuous part S( T ) defined on a semi- 
discrete curve in the space of the discrete independent variable N and the the continuous 
variable r. This is most easily visualized in the case where we have a piecewise linear curve 
as in Fig. 01 where in the horizontal direction we have step-wise jumps over units, while 
vertically the curve is continuous. The closure relation guarantees, nevertheless, that any 
such semi-discrete curve may be "locally" deformed without changing the action functional 
consisting of sums of the form (|5.1|) over the horizontal iterates and of the form (|5.6p in the 
vertical segments . 

T 



> N 



Figure 4: The semi-continuous curve. 



Variational principle for a semi-continuous 1-form structure: In (|5.7p . we obtain 
the semi-continuous closure relation. It is interesting to see how the variational principle 
works in this case. We will end this section with a scheme for deriving the closure relation. 
Now, suppose we have a curve T living on the semi-continuous (N — reconfiguration shown 
in Fig. ([5]) . The action associated with the curve T can be written in the form 

S(x(N,r) ; r) = y T+5T £ (T) (x(T,N),x(T + <5r,N))dT 

+£ (n) (x(t + St, N), x(r + St, N + AN)), (5.12) 

where the first term is the Lagrangian living on the vertical (continuous) line and the sec- 
ond term corresponds to the Lagrangian living on the horizontal (discrete) line of the curve T. 

We now consider the action along the dashed curve 

S'(x(N,r) ; r') = £ (N) (x(r,N),x(r,N+AN)) 

/t+5t 
£ (t) (x(t, N + AN), x(r + St, N + AN))g?t, (5.13) 

where the first term is the Lagrangian living on the horizontal (discrete) line and the second 
term corresponds to the Lagrangian living on the vertical (continuous) line of the dashed 
curve. 
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T 



t+5t 



• x(r + 5t, N + AN) 




N+AN 



Figure 5: The semi-continuous curve T and V. 

We now impose the condition that <5S = S' — S = 0, yielding 
S§ = £ (n) (x(t,N),x(t,N + AN))-£ (n) (x(t + <5t,N),x(t + 5t,N + AN)) 

/t+<5t 
(£( t) (x(t, N + AN), x(r + (5r, N + AN)) - £ (r) (x(r, N), x(r + <5r, N))) dr. 

Using the Taylor expansion with respect to the variable r, we have 

5§ = St [£ (t) (x(t, N + AN), x(t + 5t, N + AN)) -£ (t) (x(t, N),x(t + <5t, N)) 

-J^ (n) (x(t,N),x(t,N+AN))] . (5.14) 

The term inside the bracket is required to be zero according to the condition 5§ = 0, yielding 
the semi-continuous closure relation. 



6 The full continuum limit 

In the previous section, we took the continuum limit on the discrete variable m, leading 
to a system of differential-difference equations. The full continuum limit, performed on the 
remaining discrete variable N, will lead to a coupled system of poles in the first instance, from 
which a hierarchy of ODEs can be retrieved, which is the CM hierarchy. How to perform 
this limit is inspired by the structure of the solutions of (|4.3[) . Performing the following 
computation, 

(p + fc)V«-^ = p N exp jfce + Nln j. , 



p N exp lk£ + N 



k 



k :] 



= p e p e 



p 2p 2 3p 3 

kti+k 2 t2+k 3 t 3 +k 4 t4 + ... 



V ) 

k k 2 
- + — 

p p z 



(6.1) 



we can identify the full continuum variables as the coefficients of the various powers of k, 
namely 



t t rj t rj 

tl = ? + -q + W> *2 = q - 7p, *3 = —T + r^i 

p^ p a Lp p+ 6p z 



3.2) 
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where r\ = ~. This reduces the KP eigenfunction given in (|2.4p to the following form. 

<p = ( i - iy — ^i— ) p N e -5 e /c tl+ fc^ 2+ ^3+... ; where Xi = x . _ r _ (6 3) 

The corresponding reduction to CM removes the centre of mass, which can be identified 
with the t\ flow, from the system of equations. Since = Xi(ti,t2,ts, ...) and we take the 

dXj 

ati 



assumption that ^j- 1 = 0, we find that 



dx t 1 dX.dh dX l dt 2 dX, dt 3 



dr p 2 dt\ dr dt 2 dr dt 3 dr 
1 1 dXi 1 dXi 



We also find that 



p 2 p 3 dt 2 p 4 dt 3 

1 1 dX. L 1 dXi 1 d 2 X t 



(6.4) 



p 2p 2 dt 2 3p 3 dt 3 8p 4 di 
1 d 2 X t 1 d 2 X, 



and 



6p 5 dt 2 dt 3 I8p e dt\ 

1 1 dX, 1 dXi 1 d 2 X l 



(6.5) 



P 2p 2 dt 2 3p 3 dt 3 8p 4 dt 2 
1 d 2 X t 1 d 2 X l 



6 P 5 dt 2 dt 3 18 P 6 dt 2 



(6.6) 



The full limit on the solution: Next, we would like to perform the full limit on the exact 
solution. We now rewrite (|4.3ip in the form 

Y(N, r) = (l + ±y N e i( 1+ fr\(0,0)e-i( 1+ Tr (l + ff 

-H(i + ±Y 1 -z(i + ± 



- N ( 1 +t)+f( 1 +f)" Y (0,0)e N ( 1+ t)-f( 1+ t)" 

+ + • ( 6 - 7 ) 



p j P \ P 

Using the definitions of ti, t 2 , and t 3 in (|6.2[) . we now have the solution in the form 

Y(*i, t 2l ta, ...) = e -A(ti-O+A 3 i 2 -A 3 t3 +... Y ( 0) Qj )e A( tl -«)-A 2 t2 +A 3 t 3 + ... 

-(t!-0-2A< 2 -3A 2 t 3 + .... , (6.8) 

which is a function of time variables (t\, t 2 , t 3 , ..). The positions of the particles Xi{t\, t 2 , t 3 , ...) 
can be computed by looking for the eigenvalues of (|6.8p . Since, the term t\ — £ represents 
the centre of mass flow, the explicit expression of the solution for the CM can be obtained 
from the secular problem for the matrix 

X(0, 0) - 21.(0, 0)f 2 - 3i 2 (0, 0)i 3) (6.9) 

where X(0,0) = U _1 (0, 0)Y(0, 0)17(0, 0) and £(0,0) = 17 -1 (0,0)A17(0,Q). We note that 
(I6.9[) is identical to what was obtained in [5] up to the time variable t 2 (involving a scaling 
of the variable). The solution (|6.8[) involves iV-time flows for the CM system. The next 
solutions in the hierarchy can be generated by pushing further on with the expansion. 
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The full limit on the equations of motion: We now would like to see what would re- 
sult from taking the limit on the equations of motion. We know that the combination of eqs. 
(|4.21[) and (|4.25|) gives us the equations of motion in the skew limit case. We look at the full 
limit for each of these. 

Using eqs. ({0)1 . (jo3|) and we find that eq. (TOTj) gives 

/ 1 \ 1 / dXi \ 2 lSXi A 1 , . 

°ta) : z(^r) +3^^(x-^F^ ' ( 6 - 10a ) 
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(X-X,) 2 (X*-X 



, (6.10b) 



1\ 1 fdXiY 1 9 2 I, 3 fdXi \ dX t ldX,d 2 X. 



J I i \ »f„ / 3 9i 2 ^3 i\dt 2 J dt 3 8 9t 2 at 2 

i /ax\ 4 3axax 8 A / 3 
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a< 2 ; 4 at 2 du ^\--\, .v, 

dXj_ o dXj dXj \ 

dt 3 z gfa 9*2 
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(X-X,-) 2 (X-X,) s 



= , (6.10c) 



1 \ 1 /ax a 2 lax x 



at 2 ; 3 at 3 ^(.v, v , . 

i\ ia 2 x 4 2aXiax, 1 /ax\ 3 iax 4 
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8 at 2 dt\ 



1 /axA 4 3 ax ax 



dXj 9 dXj_ dXi 

du . Z at 2 " 1 " 9*2 



(X-X,-) 2 (X-Xj)^ 
The difference between (|6.10[) and (16. lip gives the following 
♦The leading term of order 0(l/p 3 ) gives us 

a 2 x, : 



= 0. (6.11c) 



32- N 



dtl fZ&i-Xj) 3 ' 
which is the equation of the motion for the continuous CM system 
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♦The term of order 0(l/p 4 ) gives us 



3dt 2 dt 3 + 4 dt 2 9t?, (Xi-X,-) 3 ' 1 J 



We now use eq. (|6.12p to simplify eq. f|6. 13f) . and we get 



„ 2 I? ^ , Ml 

This equation represents the next equation of motion beyond the usual continuous CM in 
the hierarchy. We will stop at this equation, but we can actually get the higher terms of 
the equation in which the variable £4 and higher order time- flows must be taken into account. 

Note finally that summing (|6.10p and (|6.1ip gives the following equations 

1 \ 1 / dXi \ 2 ldX, ^ 1 



3= 

3V* 



/ 1 \ 2dx % dXi 1 /aiA 3 1 dXi A g| 







1 / ax, \ 3 3 / dX; \ 2 dx, 1 / dX; \ 4 3 ax< ax,- 



p 4 J ' 3\dt 3 J A\dt 2 ) dt 3 2\dt 2 J 4 <9i 2 9*. 
3 gtg 

(X - X,-) 4 (X - X,-) 



N f dX 3 \ 



These equations are the constraints for the full limit case which will play an important 
role in the next section. The occurrence of these constraints seems to be related to the 
superintegrability of the CM system, cf. [5], and their validity can be checked on solutions. 



7 The continuous Lagrangian 1-form and its closure re- 
lation 

In the previous section, we obtained directly the equations of motion of the continuous CM 
from the full limit. In this section, we obtain the continuous Lagrangians from the full limit 
as well. We will restrict ourselves to the case of 3-particles involving time flows t% and t 3 , 
to establish the simplest example of a Lagrangian 1-form for the CM system. Both the 
Lagrangians &( T ) an d &(N) produce the hierarchy of the continuous CM Lagrangians corre- 
sponding to the equations of motion in the previous section. Rather than performing the 
full limit on the individual Lagrangians, we prefer to take the limit on the action itself. 

The full limit on the action: In Section 5, we performed the skew limit on the dis- 
crete action. In the present section, we will continue to perform the full continuum limit on 
the semi-continuous action. We now take the action to be of the form 

S[x(N,T);r] = f 2 drL iT) (x(N,T),x(N 1 T))+Y,^m(x(N,T),x(N + l,T)) , (7.1) 

Jti n 

where the first term belongs to the vertical part and the second term belongs to the horizontal 
part of the curve given in Fig. (jU). 
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Using anti- Taylor expansion, the action now becomes 



I-T2 rm 

§[x{N,r);T} = dTL (T) (x(N,r),x(N,T))+p dr)L m {x(N,T),x(N + l,r)) , (7.2) 

where we do not need to take into account the boundary terms coming from the expansion, 
because they are constant at the end points and do not contribute to the variational process. 
We now perform a change of variables (t, rj) i— > (t 2 ,t 3 ) by using 



dr = 2p 3 dt 2 + 3p 4 dt 3 , (7.3a) 
drj ~ —6pdt 2 — 6p 2 dt 3 , (7.3b) 



and also expand the Lagrangians with respect to variable p. We obtain 

rU[2 \u r („n , ^ dx(t 2 ,t 3 ) 



SI , . v Fl F W A+ , ( . , dx(t 2 ,t 3 ) \ 

S[a!(*2,t3);r] = / dt 2 L( t2 ) [x(t 2 ,t 3 ), 

Jt 2 {i) V ot 2 j 

r („n * \ 9x(t 2 ,t 3 ) dx(t 2 ,t 3 y 

*3(1) 



r ( . x dx(t 2 ,t 3 ) dx(t 2 ,t 3 ) \ 
dt 3 L {t3) ^ x (t 2 ,t 3 ),—^— > —^—j , (7.4) 



where £(t 2 ) an d £(t 3 ) are given by 



ifdx^- x 



of which the Euler-Lagrange equation 



gives exactly eq. (|6.12p and 



= , (7.6) 



We see that the Lagrangian £(t 3 ) contains derivatives with respect to two time flows t 2 and 
t 3 . We observe that the equation of motion (|6.14[) requires the Euler-Lagrange equation in 
the form 



ac (ts0 a / ac (t3 



5X 4 dt 3 \d{^) t 

Furthermore, we find that 



. (7.8) 



2 N 



which is identical to eq. (|6.15a[) . 

Here we obtained the hierarchy of Lagrangians for the CM-model through the full contin- 
uum limit. Obviously, higher Lagrangians in the family can be generated by pushing further 
on with the expansion. 

The full limit on the closure relation: Using the results regarding the limits of the 
Lagrangians £>M an d £(n) m the previous section, and also the relations of eq. f|6 . 2[) . we can 
write 

dL W , , 1 d£ (t2 , t3 ,...) 1 ^( t2 ,t 3 ,...) 

~d^^~p^~dt 2 + p^~dt 3 ■ (7 - 10) 
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Substituting eq. (|7. 10[) into the closure relation eq. (|5.7I) together with the full limit of 
Z( T ) — £(t)j we nn( i that the leading term is of order 0(l/p 3 ), and is 

This equation represents the closure relation for the continuous case of the Lagrangian 1-form. 
We are now taking into account not only the usual CM Lagrangian, but also a higher-order 
Lagrangian; these have up until now not been considered. The verification of the closure 
relation (|7.11[) is presented in Appendix [D] 



8 The variational principle for forms 



In the previous section, we obtained the Lagrangian hierarchy of the CM system along with 
the corresponding Euler-Lagrange equations. In this section, we are interested in investigat- 
ing the variational principle for continuous Lagrangian 1-forms. Here we limit ourselves to 
the first two time flows in the hierarchy, t 2 and £3, corresponding to a 3-particle system in 
order to exhibit the simplest example. Let us define the time variables to be functions of a 
new variable s, so that fj = U(s), associated to the curve T : (t 2 (s), t 3 (s)) as in Fig. 
We now define the Lagrangians to be 



(*2) 



(* 3 ) 



-C(ij)(x(i 2 , t 3 ), x t2 (i 2 , t 3 ), x t3 (t 2 , t 3 )) , 
<C(t 3 ) (x(*2, t 3 ),x t2 (i 2 , t 3 ), x t3 (t 2 , * 3 )) , 



where 



<9x c'x 

*t 2 (t 2 ,t 3 ) = — and x t3 (i 2 ,t 3 ) = ~— 



<9x 



0t 



a< 3 



.la) 
.lb) 



The action evaluated on the curve T can be defined as 



S(x(t 2 ,t 3 );r) 



(£( f2 )dii + L^dt 2 



^ (t2) (i 2 ( s ),t3( s ))^+£ (t3) (i 2 ( s ),t 3 ( s ))^^d S . (8.2) 



i 3 




•*2(so),t3(sO 

(a) 



*2 



i 3 



*2(si),*3U 




*2 



Figure 6: The deformation of the continuous curve T . 

We now perform variation on the curve of time variables t 2 and t 3 in the following way. 
Suppose we have an action § defined on a curve F, and the curve is deformed t 2 ^ ^2 + <5^2 
and t 3 1 y t 3 ~\~ St 3 in Fig. (JS}}), keeping the end points fixed, so we have the conditions 

St 2 (s ) = 8t 2 ( Sl ) = 0, 5t 3 (s ) = 5t 3 ( Sl ) = 0. (8.3) 
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The action §' is denned as 

§' = §(x(t 2 + 5t 2 ,t 3 + st 3 )-r r ) 

Extremizing the action, we have <5S = 0, where SS is defined by 

S§ = S(x(ta + 6t 2 , t 3 + St 3 );T') - S(x(t 2 , t 3 ); r). 
Using the Taylor expansion, we obtain 



(8.4) 
(8.5) 
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Using the fact that the end points are fixed, the first line of eq. 
use 

dL/ 



J (*2) 



ds 



dL 



its) 



to obtain 
5§ 







( ds 


st 2 ( 


J s 





ds 



d£(t 2 ) dt 2 <9£(t 2 ) dt 3 

dt 2 ds dt 3 ds ' 
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(8.6) 

is zero. We can also 

(8.7a) 
(8.7b) 
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Here we must have two copies of the closure relation 

&C(t 3 ) _ 9£(t 2 ) 



<9z 2 



az. 



(8.9) 



in order that 5§ = 0. The closure relation (|8.9[) guarantees that any such continuous curve 
may be locally deformed without changing the value of the action functional corresponding 
to the time variables t 2 and t 3 . 

Next, we would like to perform a variation on the field variables x. We now fix the curve 
r on the t 2 — 1 3 plane and examining the corresponding curve £r of the variables x living in 
the x — t 2 — t 3 space as in Fig. ((7^,) . The action can be expressed in the form 

s(x,y 2 ,y 3 ; £ r) 

£(t 3 ) (x(i 2 , Z 3 ), x t2 (Z 2 , t 3 ), x t3 (t 2 , t 3 ))dt 3 + £ (t2 ) (x(Z 2 , t 3 ), x f2 (t 2 ,t 3 ),x t3 (t 2 ,t 3 )) dt 2 ) 

dt 3 



ds [ L (t3) (x(t 2 (s),t 3 (s)),x t2 (t 2 (s),t 3 (s)),x t3 (t 2 (s),t 3 (s))) — 



+ £ (t 2 ) ( x ( t 2(s),t 3 {s)),x t2 (t 2 (s),t 3 (s)),x t3 (t 2 (s),t 3 (s)) 
We now introduce 
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ds 
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dU 



Sx t2 — + Sx t3 — 



dt 2 dt 3 
5x t —-5x t3 — 



(8.10) 

(8.11a) 
(8.11b) 
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Figure 7: The curves T and £r in the x — t configuration. 



where y is a new independent variable, and the combination of these two equations yields 

(8.12a) 
(8.12b) 
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Performing the variation on the action, we have 
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Using dHH21), (jglU)) becomes 



si 

ds 

so 

1 d 



dL( to ) dt 2 9L( U \ dt. 



(*3) 



<9x ds 



dx 



<5x 



dL 



C*2) 



<9£ 



(*3) 



dt 2 /dsdL( t2) dt 3 /dsdL {t3 



dx t2 



dx t . 3 dt 3 /ds dx t3 dt 2 /ds dx t2 



9£ (t2 ) 9£ (t3) dt 2 /dsdL {t2) dt 3 /dsdL [t3) 



dx t2 chc t3 dt 3 /ds dx t3 dt 2 /ds dx t2 
Integrating by parts the second line in (18.14[) with the boundary conditions 

Sx(*2(so)) = <yx(t 3 (*o)) = , «x(*a(«i)) - <M* 3 (si)) = , 



we obtain 
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The condition 5S 
c)L 



is satisfied once 
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= 0. 



(8.16) 
(8.17) 



We consider eq. (|8.16l) to be a generalized Euler-Lagrange equation containing two separate 
Euler-Lagrange equations, and the equation (|8.17[) is the constraint resulting from the vari- 
ational principle. 




— """-' •(t2(gl),*3(gl)) 



Figure 8: The deformation of the curve T. 

To summarize, performing the variation on the curve T with respect to the time variables 
on the t 2 — 1% plane, the computation leads to the closure relation (|8.9p . By fixing the curve 
r and varying the curve £r, the computation leads to the Euler-Lagrange equations corre- 
sponding to the Lagrangian 1-form. However, the closure relation allows the deformation of 
the curve T to a simpler curve, as in Fig. ([S]). If we now are on the curve T 2 where the time 
variable t 2 is a constant leading to — 0, we obtain 



ac 3 9 

dx dt 3 
dL 3 



^3 



dx t . 



= 0. 



(8.18a) 
(8.18b) 



If we are on the curve Ti, where the time variable is a constant leading to ^ = 0, we 
obtain exactly the same equation (|8.18p with interchange of indices 2o3. This is entirely 
consistent with the results obtained in the previous section. 

The Euler-Lagrange equations (|8.16[) . the constraint equations (18.17[) and the closure 
relation (I8.9[) together form a consistent set of equations, representing the multidimensional 
consistency of the system. 



9 Concluding remarks 

In this paper, we present what we believe to be the treatment of the Lagrangian structure 
for the Calogero-Moser system which best captures the essential integrability characteris- 
tics. It is well known that the continuous CM system is Liouville integrable [2] and also 
superintegrable [5] and this remains the case on the discrete-time level [351 US]- However, 
in this paper, we focused on another aspect of the integrability, namely multidimensional 
consistency, meaning in this context of finite-dimensional systems the existence of commut- 
ing flows. We believe that in the present paper, we have given convincing arguments for 
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the assertion that the proper Lagrangian structure should be the one in terms of Lagrange 
1-forms. Intriguingly, this is most manifest on the discrete-time level where it follows directly 
from the relevant Lax equations. In the continuous case, the relevant Lagrangians (apart 
from the obvious Lagrangian for the second order CM flow) are more difficult to establish, 
but they follow through systematic continuum limits performed on the discrete Lagrangian. 
However, the latter is quite a subtle computation, and it is here that the connection between 
the CM system and the KP system (the former arising from the pole reduction of the lat- 
ter) is essential. The construction in this paper starts from the semi-discrete KP equation, 
and it is the dependence of this equation on (lattice) parameters which guides the choice 
of continuum limits. We believe that the CM system forms a first important example for 
the study of the new variational principle which applies to the case of Lagrangian 1-forms, 
cf . [TH [2TJ [221 [31] , as preparation for similar structures in the case of higher Lagrangian 
multiforms. 



We finish with a few observations. 

• Since the construction is based on the connection with a semi-discrete KP equation, 
it would have been more satisfactory if a Lagrangian for that equation was at our 
disposal. However, at this stage, such a Lagrangian seems to be elusive, even though 
it exists for the bilinear discrete KP case |22) . 

• Another issue is the connection with the Hamiltonian hierarchy associated with the 
Calogero-Moser system. Usually, the connection between the Hamiltonian and the 
Lagrangian is given through the Legendre transformation, which for the second flow 
in the CM hierarchy is the standard one. However, when we move to the higher 
Hamiltonians in the hierarchy, these are no longer of Newtonian type, and hence lead 
to complicated expressions when implementing the Legendre transformation. What 
the treatment in of Section 7 reveals , however, is that the higher flows fit naturally 
in the 1-form structure through mixed, but polynomial, Lagrangians in terms of the 
higher-time derivative. 

• The natural framework for quantization of the discrete-time CM model seems to be 
that of the path integral, which can be interpreted as a quantum version of the least- 
action-principle. One may conjecture that the quantization scheme exploiting the new 
Lagrange structure could lead to a rigorous approach to constructing the relevant path 
integral. 

Our motivation in this paper was to present the most simple example of a Lagrange 1-form 
structure, hence our restriction to the rational case. However, most of our results can be 
readily expanded to the trigonometric and hyperbolic and also elliptic cases. The formulae 
for those cases are presented in Appendix [B] A natural question is how our results extend to 
the relativistic case, i.e., the Ruijsenaars-Schneider (RS) model. We intend to address that 
problem in a future paper. 



A The construction of the exact solution 

In this appendix we review the construction of the exact solution (|2.12[) . which follows a 
procedure similar to the continuous case, cf. e.g. [S]. The basic relations following from the 
Lax pair, i.e., (|2.1ip together with the definitions (|2.7p and (|2.8[) . lead to: 

(L-M)E = , (A.la) 

E{L-M) = 0, (A.lb) 

XM - MX = E , (A.lc) 

XL LX = I-E, (A.ld) 

where X = X(n,m) = 53i=i x i( n : m )Eu is the diagonal matrix of particle positions. We 
concentrate first on the part of the Lax pairs associated with the dynamics in terms of the 
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variable n. From the Lax equation (|2.11[) . we may write 

M = UU~\ and L = UAlU' 1 , (A.2) 
where U = U(n, m) is the matrix used to diagonalize L. Then (|A. 1|) becomes 

U~ X E = A^U^E, (A.3a) 
EU = E UA L , (A.3b) 

Y-Y = -JJ~ X EU , (A.3c) 
YA L -A L Y = U-\I-E)U, (A.3d) 

where Y = U~ 1 XU, which leads to the following expression 

YA L - A L Y = I . (A.4) 

Noting the invariance of A^ under the discrete-time shift, we can easily solve (|A.4|1 . and its 
general solution is 

Y(n, m) = Al n Y(0, rn)A n L - nA^ 1 , (A.5) 

with Y(0,m) determined from the initial data X(0,m). A similar analysis can be applied 
to create the solution associated with the " ^ " shift. 

Conversely, we can start from a given N x N diagonal matrix A with distinct entries, 
and an initial value matrix Y(0, 0) subject to the condition that 

[Y(0, 0) , A] = J" + rank 1 , (A.6) 

where [, ] denotes the matrix commutator bracket. Let J7 -1 (0,0) be the matrix that diago- 
nalizes Y(0, 0), i.e., such that 

F(0,0) = (7 _1 (0,0)X(0,0)(7(0,0)) , X(0, 0) = diag(zi(0, 0), . . . , x N (0, 0)) . (A.7) 

If the eigenvalues of 5^(0, 0) are distinct (which we can take as an assumption on the initial 
condition) then C _1 (0, 0) is determined up to multiplication from the right by a diagonal 
matrix times a permutation matrix of the columns. (Fixing an ordering of the eigenvalues 
Xj(0,0), f/ _1 (0, 0) unique only up to multiplication by a diagonal matrix from the right). 
We can fix C7 _1 (0, 0) up to an overall multiplicative factor by demanding that 

[Y(0, 0) , A] = I — C/ _1 (0, 0) E U(0, 0) . (A.8) 



Next we consider the matrix function given by 

Y(n, m) = A- L n A K m Y(0, 0)A£A£ - nA£ J - mA^ 1 , (A.9) 

with Al = pi + A and A#- = ql + A. Let C7 _1 (n, m) be the matrix diagonalizing Y(n, m); 
by an appropriate choice of an overall factor (as a function of n and m) this matrix can be 
fixed such that it obeys: 

U- 1 (n,m)E = Al n A K m U- 1 {0,0)E , E U(n, m) = E 17(0, 0)A£A£ , (A.10) 

and 

[Y{n, m) , A] = I — C/ _1 (n, m)E U(n, m) . (A.ll) 
From the expression (|A.9[) we can now derive the relations 

A L Y Y A L = I , (A.12a) 
A K Y - YA K = I , (A.12b) 
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with the usual notation for the shifts in n and m over one unit. Together with the relation 
(jA.lip this subsequently yields: 

Y Y = TJ EU , Y-Y = -U~ 1 EU. (A.13) 

Reversing these relations by rewriting them in terms of X(n, to) = U(n, m) Y(n, to) U" 1 (n, to) 
and now defining the Lax matrices by 

L:=Uh L U- 1 , K:=UAU\ (A.14) 

together with 

M:=UU^ , N-^UIJ- 1 , (A.15) 

we recover the relations: 

[X, L] = [X, K] = I-E , XM MX = XN- NX = E , (A.16) 

which determine the matrices M and N as functions of the Xi(n,m) as well as the off- 
diagonal parts of the matrices L and K. From the definitions of L and K we have that 

L — K = (p — q)I . 

Furthermore, from (|A.12I) we obtain 

L X M M X L M , KXN-NXK = -N, (A.17) 
which, when combined with the latter relations of (|A.16[) . yield 

(lM-ML)(X-E)=0 , (KN- NK^ (X-E) = . (A.18) 
On the other hand, using the first of the relations (|A.16[) we also obtain 

X (1m - ML} = E(L- M) and X (^K N - N K^j = E (K - N) , (A.19) 
as well as 

f LM -ML)X = {L-M)E and (kN-Nk)x = (K-N)E. (A.20) 



From the relations (|A.18|I , (|A.19|I and ([A. 201) it follows that the Lax equations hold and their 
form is determined up to the diagonal part of the matrices L and K . 

B The elliptic and trigonometric discrete-time Calogero- 
Moser 

In this section, we show that there is a connection between the time-part Lax matrix and 
the Lagrangian for the trigonometric/hyperbolic and elliptic cases of the discrete-time CM 
system, similar to that which we have established for the rational case. 

B.l The elliptic case 

The Lax matrices in this case read [23] 

N I N N \ N 

L ^ = E E c(* - ?i) - E c(?< - qj) ) E u - E - «>-)3y - (B-i) 

i=l \j = l j^i J i^j 



and 



N 

M eU = J2 'i'. ''/- '/ ; , (B.2) 
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where ((x) is the Weierstrass zeta function and & K (x) — ^^t^t, where a(x) is the Weier- 



t(k)<t{x) 

strass sigma function. The Lax equation LM = ML produces the equation of motion 



N 



N 



3=1 V 



(B.3) 



and the corresponding Lagrangian is 



N 



N 



£>ell =~ log|cr(g; -q~j)\ +^log|cr(9i - qj)\ 

j,J=l if 3 



(B.4) 



The determinant of the matrix M e ui p ti C is 
det|M e «| = $ K (£)(7(£) 



el - <7i) 



(B.5) 



where S = 53i=ite ~ an d then we can also write 

AT 

log|det(M eii )| = log |$„(E)o-(E)| + $^(log |cr( 9i - qj)\ + lo S k(S 

AT 

- ^log|<7@- gj )| . (B.6) 

ij 

The action of the system can be expressed by considering the chain product of the matrix 



Jell 



log 



dct Yl M eii{n) 



\n—+oo 



N 



N 



Y log|$ K (Ey(E)|+^log|c7fe-g,)| -^logla®-^ 



>f-j 



^(£ e «+l0g|$ K (E)<7(E)|) 



(B.7) 



The last term in eq. (IB. 71) is related to the motion of the centre of mass of the system and 
it can be separated from the relative motion. 

B.2 The trigonometric case 

The Lax matrices in this case read [U] 



N I N 



N 



N 



'"■■■> = g yC m " i::,/ ' !/ ' ' § C ° th( * - * } J ^ - g iinh , mu!„ k . 



sinh(q; - ^ + . 

(B. 



and 



N 



M 



trig 



ESI 
si 



sinh(cfj - 9j + k)^. 



^-^ sinh(gi — qj) sinh(«) 
The Lax equation LiVf = ML produces the equation of motion 



(B.9) 



A? 



A r 



^ cothfe + coth(% - £j) - 2 ^ coth(% - qj) = , (B.10) 

3=1 V / ]fi 
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and the corresponding Lagrangian is 



N N 
&trig = - 2J l0 S l silm (<7i ~ 9j)\ + X] l0g l sinh (* _ Qj)\ 

The determinant of the matrix Mtrig is 

sinh(£ + k) IL<, sinh(g. ( - qj) sinh(<fc - q~j) 



det \M 



trig 



and then we can also write 



log det (M trig )\ = log 



sinh(/t) 



sinh(£ + k) 



n« sinh® - qj) 



(B.ll) 



(B.12) 



JY 



sinh(K) 
^log|sinh(gi - qj 



i<3 



N 



5^ (log |sinh(% - + log |sinh(ft - q 3 )\) 

(B.13) 



The action of the system can be expressed by considering the chain product of the matrix 

Mtrig 



^trig 



log 



det ( J] M trig (n) 



= ( &trig 



sinh(E + k) 



s'mh(n) 

sinh(E + k) 
sinh(fi;) 



N N 

log Isinfe - qj)\ - ^ log \cr{q, - ^) 



(B.14) 



Again the last term in eq. (|B.14[) is related to the center of mass motion of the system, and 
it is not going to effect the equations of the relative motion. 

C The derivation of the Lagrangian Lr T \ 

In this section, we show the derivation of the Lagrangian Li T \ (15.8[) from the skew limit of 
the Lagrangian £( m ) (I3.10b[) . We recall the Lagrangian Jdi m \ again here 



£ (m) =log|det(JV)| +g(H-S) 



(C.l) 



Performing the skew limit, we obtain 



(m) 



log I dct(M - eA)| + (p - e)(S - 3 - eE) 

log I det (M) I + log I det(l - eM ^A) + ... + p(E - 3) - e(S - 3) - peE 

... - e(S - 5) -peE 



£ (N) +log|l-eTr^M _1 A 

£(N) + e£( r ) + 



where 



r) = -Ti^M^a) -E + E-pE 



The inverse of the matrix M is given by 



N 



m- x = -E 



E- ■ 



(C.2) 
(C.3) 

(C.4) 



ij=l 
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where 



We now can show that 



*(*) = "|\ (Xt * } , (C5a) 
Ui^i -xi) 

= "If 3 "^ - (C5b) 



Tr(M-'S)-f»^. (C. 6 , 

1,1 = 1 J 



We now consider the identity 

which wc differentiate with respect to £ at a given vafue £ = x,- , yieiding 



We now differentiate (|C.8j) with respect to Xj, yieiding 



Using (|C.9p . we can rewrite (|C6[) in the form 



3 1 -eS. 



a i 



Tr ( M _1 A ) = Vl*,)^ = 



N 



then the Lagrangian &< T ) can be expressed in the form 



E^tI- 111 ( c - 10 ) 



N ~ N ~ 

Xi x X 



. . Xj X.y . . Xj X 



D The direct proof of the closure relation ( 17.111 ) 

In this section, we show that the closure relation (17. lip holds on the equations of motion. 
Here we recall the equations of motion for the time variables ti and £3 



N 



d 2 x t 

dt 2 dt 3 b 2_. (x ._ x .)3- {U - Z) 
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The Lagrangians corresponding to (|D.1|) and (|D.2I) take the form 

<*> - tlM'f^m 1 ^' (D - 3) 



respectively. 
We find that 



3r iV o(dX t \2 odXi 



and 



We see that the first terms of eq. (|D.5|) and eq. (|D.6I) are identical. The remaining work is 
to show that the last term in eq. (|D.6I) is zero. Using eq. (ID.1[) . we may rewrite the last 



term in eq. (ID. 6ft with the help of the identities: 

AT d 2 Xj N N 

1 dt 2 = \ " \ " 

8 Zv (Y; _ V.-V2 Z^ Z^ 



AT j JV 1 

We see that the first term is an antisymmetric function and hence vanishes, while the second 
term 

N 1 

= j2 f 1 i \ 2 i i 



N 

E 



\Xi — Xj Xi — X k J (Xj — Xk) 2 Xi — X k 

1 



(x i -x j ) 2 (x j -x k ) 2 (x i -x k ) 

1 2 



(Xj - x k ) 2 (x t - x k f (Xi - x 3 )(x % - x k y)(x 3 - x k 



(D. 



The first and third terms are the antisymmetric functions, hence vanish leaving the middle 
term which is actually the opposite of the term on the left hand side, i.e., 

N 1 N 1 

. E fe (x i -x j y(x i -x k Y = . (x k -x 3 y(x k -x t y = ■ (D ' 9) 
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